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Abstract 

Improved results for Bk are discussed. Scaling corrections are argued to be 
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1 INTRODUCTION 

I discuss here the improvements made by the "staggered" group (Gupta, Kilcup, Patel and 
myself) since LATTICE 91 in the calculation of 

' l{K\s^,^,d\0){0\s^,^,d\K) • ^ ^ 

In 1991 I quoted results using lattices at /5 = 6, 6.2 and 6.4 The major uncertainties 
were: 



• Linear vs. quadratic extrapolation in the lattice spacing a. These two choices led to 
Bk = Bku-^'^ = 0.66(6) and 0.78(3) respectively, with the errors being statistical. 

• Perturbative corrections were not included, because of the uncertainty in the choice of 
9- 

• Contamination from excited states. 



• Quenching, although there was preliminary evidence that this was not the dominant 
source of error |]. 

• Extrapolation from degenerate quarks (mi = 1712 ~ ms/2) to the physical kaon (mi = 
rris, 1712 ~ 0). From the quenched data, we estimated that this extrapolation increased 
Bk by 3%, an estimate included in our results for B^- The error in this estimate was, 
however, large. 

A number of these uncertainties have now been substantially reduced, both due to our 
work and that of the Kyoto-Tsukuba group. 



2 SCALING VIOLATIONS ARE 0[a^) 

According to the standard lore, the staggered fermion action is good to O(a^) (up to loga- 
rithms, which will always be kept implicit). I sketched the perturbative argument for this in 
Ref. |1|, and presented some supporting numerical evidence from the spectrum. By contrast, 
for matrix elements of external operators the corrections are expected to be of 0(a). For 
example, the operators appearing in both numerator and denominator of B^ (Eq. |I]) have 
0(a) terms in their tree- level perturbative matrix elements f^. These 0(a) terms turn out 
to have the wrong flavor to contribute to the matrix element in Bk-, but in 1991 it seemed 
possible that terms of 0{g'^"'a), from n-loop diagrams, might contribute. 

In the following I sketch an argument that they do not the scaling corrections to the 
matrix elements in B^ are of O(a^) to all orders in perturbation theory. The argument is 
an application of Symanzik's perturbative improvement program 0, 

It is useful to begin by demonstrating that the staggered action is already "improved", 
i.e. has no corrections of 0(a). In the notation of, e.g., Ref. 0, |^, the action is 

s = T. xiy) [il,^I)D, + iTWl)m]xiy) , (2) 
y 
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and is invariant under translations, rotations, spatial inversions and charge conjugation. 
When m — 0, it is also invariant under the axial transformations 

X exp[za(75 ® i^)]x , X ^ Xexp[ia(75 (g) ^5)] . 

In order to improve the action, one adds operators of d = 5, with coefficients adjusted 
order by order in perturbation theory so as to cancel 0{a) terms in correlation functions. 
(In fact, except in scalar theories, only on-shell quantities are improved.) That this improves 
all on-shell quantities at once is the non-trivial assumption, shown by Symanzik for scalar 
theories. The d = 5 operators must have the same symmetries as those in the original action. 
Ignoring axial symmetry, the allowed operators are |^ 

m^xil ® I)X , X{I ® I)Dlx , mxilf, ® I)D^x , xkly^v ® l)F^,x ■ (3) 

However, none of these operators is consistent with the axial symmetry — treating m as a 
spurion field one can show that the bilinear must either contain an even number of links, and 
be multiplied by an odd function of m, or contain an odd number of links and be multiplied 
by an even function of m. Since there are no operators available to improve the action, it 
must already be good to O(a^). 

Now I proceed to Bk- I will discuss the matrix element of the four-fermion operator Ob 
in the numerator of Eq. |l]; a similar argument works for the simpler matrix elements in the 
denominator. To improve a matrix element one must not only improve the action, but also 
improve the operator itself. Since the staggered action is already improved, one needs only 
consider the operator. I assume, following the second paper in Ref. 0, that improvement 
can be accomplished for all on-shell matrix elements by adding = 6 and 7 operators to 
the original operator, with coefficients determined order by order in perturbation theory. 
(Mixing with lower dimension operators is forbidden by the flavor structure.) I assume 
further that these operators must have the same symmetries as the original operators. 

Let Ocont tie a vector of continuum = 6 operators. At tree-level these match onto 
a corresponding vector of lattice operators 0\g^. These lattice operators will mix with all 
others of (i = 6 having the same symmetries, so one must extend the vector 0\^^ to include 
all the possibilities. I assume that the continuum vector is similarly extended. Let 0\g^ be 
a vector containing all operators with c? = 7 at tree-level having the same symmetries as the 
0\^^. The assumption of all-orders improvement is then 

+ 0(a^)) = c(/)OL + d{g^)aOl, , (4) 

in the sense that lattice matrix elements of the operator on the r.h.s. equal those of the 
continuum operators up to corrections of 0[c?\ The coefficients c and d are matrices, c is 
square, and is the identity at tree- level, d is rectangular, and begins at 0(1) in general. 

We do not want to have to calculate c and d to all orders. Thus we want an expression 
for 0\^^ in terms of continuum matrix elements alone. This requires relating 0\^^ to the 
corresponding continuum operators 0\^^^. The assumed form is 

aOLt(l + 0(a)) = ~c{g')Ol, + d{g')aOl, , (5) 

where d is a square matrix which is the identity at tree- level, c = Olg"^) represents the fact 
that d = 7 operators mix back into d = 6 operators. 
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Reorganizing these equations we find 

Ol, = (1 + 0(a^))(c - dd-'cr' (oLt - dd-'aOl^,) . (6) 

This shows that the 0(a) terms in matrix elements of Of^^ can be obtained from the con- 
tinuum matrix elements of Oconi- The perturbative matrix on the r.h.s. is of the form 
1 + Oi^g"^) — multiplying by its inverse calculated to, say, 1-loop, removes 0{g^) corrections 
to the matching between continuum and lattice matrix elements. 

I now apply this equation to Ob- One must transcribe the continuum operator onto the 
lattice, and then write down the (long) hst of operators Of^^ and 0]^^. Various choices of 
lattice operator, all agreeing at tree-level, have been used. We use Landau-gauge operators 
without gauge links on 2"^ hypercubes, and smeared Landau-gauge operators on 4^ hyper- 
cubes (these have d = Oi^g^)). Gauge-invariant operators have been used in Ref. P]. 
The argument is the same for all these choices, because they behave the same way under 
the relevant symmetries: the hypercubic cube, and the separate axial rotations of the four 
fermion fields. 

The crucial point is this. With staggered fermions the continuum theory has four degen- 
erate versions of each quark, and a corresponding flavor symmetry. The continuum operators 
of interest have flavor ^5 x ^5. It turns out, however, that none of the d = 7 operators has 
this flavor ||^. Thus, if we take the matrix elements between a K and K both of flavor ^5, 
the contributions of Oconi vanish identically. Thus there are no 0(a) corrections to these 
particular matrix elements: they are automatically improved. 

There are other operators in O^^^^^ having flavor ^5 x ^5, but these are multiplied by 
coefficients of 0((?^) {0{g'^) if one uses one-loop matching). 

A concern with this argument is that c and c might contain non-perturbative parts of 
0(a). This will not matter, however, as long as the symmetry properties are retained. 



3 OTHER IMPROVEMENTS 

We use the same set of lattices as in Ref. but we now have results with two sets of 
operators: Landau-gauge unsmeared and smeared. Both have O(a^) scaling corrections, but 
they should agree in the continuum limit. This provides a consistency check. 

We have now included one-loop perturbative corrections. Patel and I have calculated 
these for both the original and smeared operators]^, |^, the former results being in agreement 
with those of Ref. HSIl. The results are of the form 
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To extrapolate to the continuum, we choose a fixed scale: /i = 2 GeV. 

Lepage and Mackenzie have shown that perturbative corrections are reliable if one uses 
the correct expansion parameter []TU| . For the coefficient of SO, we use g^ determined from 
Tt{U) in Landau gauge, which yields g^j = 1.82,1.66,1.54 for /? = 6,6.2,6.4. For the 
coefficient of logarithm, which represents the effect of loop momenta between n/a and fi, we 
use either gfj, or the value obtained by running from gfj = 1.82 at /? = 6 to = 2GeV using 
the 2-loop Nf = formula, assuming that the starting scale is n/a. For 1/a = 1.9GeV at 
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Figure 1: Extrapolating Bk to a = 0. 

/9 = 6, this gives 5'^(2GeV) = 2.72. In the end we take the average of these two methods, 
and use half the difference as an estimate of the systematic error. We are in the process 
of reducing this uncertainty using the automatic scale fixing procedure of Ref. 



10 . The 



perturbative corrections are small for unsmeared operators, increasing the final result by 
~ 3%. The corrections are larger (up to 10%) for smeared operators. 

Although Bk is dimensionless, it has a weak dependence on the lattice spacing because 
of the anomalous dimension factor in Eq. |^, and through the value of the lattice kaon mass. 
We use updated values of a determined from mp. 1/a = 1.9, 2.5, 3.55 GeV for /9 = 6, 6.2, 6.4. 
Repeating the analysis using a determined from f^^ gives an estimate of the corresponding 
systematic error. 

For = 6 our bare lattice numbers are unchanged from those presented in Ref. and 
since confirmed by Ref. P]. A.i (3 = 6.2 and 6.4 our calculation is hampered by the relative 
shortness of our lattices in the time direction, which leads to contamination from more 
massive states, particularly p mesons. We use two methods of calculation each with different 
sources of contamination |11]. We now understand how to subtract these contaminations 
using the data itself |12[. To be conservative, we use the size of the subtractions as an 



estimate of the systematic error. 

An example of the extrapolation to a = is shown in Fig. [T|. The data is not good 
enough to distinguish between linear and quadratic dependence on a, although it favors 
the latter. Thus we rely on the theoretical argument given above and assume a quadratic 
dependence. An important consistency check is that smeared and unsmeared operators agree 
in the continuum limit. It turns out that they also have similar dependence on a. This is 
only true, however, after inclusion of perturbative corrections. It has been found in Ref. 
[|] that gauge invariant operators also give consistent results, including the a dependence. 
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once perturbative corrections are included. Thus the relatively large scaling violations do 
not appear to be an artifact of using gauge non-invariant operators. For our central value 
we use the average of the extrapolated results from unsmeared and smeared operators, and 
take half the difference as an estimate of a systematic error. For the statistical error we take 
the larger of the two errors. 



4 NUMERICAL RESULTS 

The preliminary result from our analysis is, in the quenched approximation and for degen- 
erate quarks, 

5i^(NDR,2GeV) = 0.616 ± 0.020(stat) 

±0.014(/) 
±0.009(sca/e) 
± 0.004(operator) 
± 0.002 {contamination) 
= 0.616 ± 0.020 ± 0.017 , 

where, in the last line, we have combined all the systematic errors in quadrature. It is 
more conventional to quote a result for the scale independent B-parameter, B^- Using the 
continuum a. evaluated at 2 GeV with = SOOMeV, and the continuum anomalous 

* MS ' 

dimension, we find 

Bk = SA'(NDR,2GeV) «,(2GeV)-''/2^ 

= 0.825 ± 0.027 ± 0.023 . (8) 

The major change from Ref. IQ is the use of quadratic extrapolation. Perturbative cor- 
rections also increase the result, by ~ 3%, and a similar increase results from the use of a 
different (and better) definition of Bk- 

using a different (and better) method of matching to the continuum Bk- 
Errors due to quenching and to the use of degenerate quarks are not included in these 
results. There are reasons to think, however, that these errors are comparable to those 
quoted above. Unquenched calculations are now possible for quark masses ~ nis/2, on 
lattices with spacing 1/a ~ 2GeV, and find results for Bk which agree within errors with 



quenched results g, |T3|. This is surprising and encouraging. It must be tested at smaller 
lattice spacings to determine whether the full and quenched a dependences are similar. 

In most quantities a more important issue would be the dependence on the light quark 
masses. For Bk, however, the dependence enters at non- leading order. If one uses the chiral 
logarithms to estimate the order of magnitude of the correction, one finds a 3% increase for 



non-degenerate quarks [14, 12 1. It is important to check this with unquenched simulations 



for rus 7^ m^. Quenched data is a not good guide because of contamination from rj' loops 
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If the result for Bk withstands further scrutiny, it will have considerable phenomenolog- 
ical impact. 
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